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In spherically symmetric, static spacetime, we show that only j = 1/2 fermions can satisfy both 
Einstein's field equation and Dirac's equation. It is also shown that neutrinos are able to have 
effective masses and cluster in the galactic halo when they are coupled to a global monopole situated 
at the galactic core. Astronomical implications of the results are discussed. 
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I. INTRODUCTION 



In Barriola and Vilenkin's global monopole solution 
|jj to Einstein's field equation, scalar fields with global 
0(3) broken symmetry are minimally coupled to grav- 
ity and the background spacetime has deficits of an- 
gle. Nucamendi et al. j| suggested that the global 
monopole solution could explain the flat rotational veloc- 
ity curves(FRVC) of stars in galaxies because its energy 
density is proportional to I/r 2 and hence it can be dark 
matter in the galatic halo. Even if some questions are 
issued H about the global monopole solution for FRVC, 
various generalized versions of the global monopole were 
studied as models for dark matter and dark energy [Q . 

Neutrino clustering was studied by some authors in- 
cluding one of us Q to explain the continuation of the 
cosmic ray spectrum beyond the GZK(Greisen, Zatsepen 
and Kuzmin) cutoff. Neutrino-antineutrino^-T/) annihi- 
lation to a Z°-boson is one possible explanation for the 
phenomena. But for v-V annihilation into a Z°-boson 
to produce super GZK air shower, a neutrino flux at 
E v > 10 21 eV with m v ~ a few eV and a significant clus- 
tering of the relic neutrino density in our galactic halo 
are required |J . 

In this article we investigate if neutrinos can coalesce 
into neutrino clouds in a curved spacetime. Consider- 
ing fermions coupled to an 0(3) triplet of scalar fields in 
the most general static metric with spherical symmetry, 
we show that only the total angular momentum j = 1/2 
fermions can satisfy both Einstein's field equation and 
Dirac's equation. When the global 0(3) symmetry of 
the Lagrangian is spontaneously broken to U(l) at the 
ground state of the scalar potential in the shape of Mex- 
ican hat, we can have a global monopole solution similar 
to that of Barriola and Vilenkin. Being coupled to a 
global monopole situated at the galactic center, neutri- 
nos are able to have effective masses and cluster in the 
galactic halo. 
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II. GLOBAL 0(3) SYMMETRIC MODEL FOR 
NEUTRINO CLUSTERING 



The action of the global 0(3) symmetric model of 
scalar fields <£> m and massive neutrinos ^„ (m, n = 
1, 2, 3) minimally coupled to gravity is given by 




S = 



(1) 
(2) 



V 

where V(Q 2 ) is a scalar potential and the last term is 
the scalar-neutrino interaction for which we assume the 
first leading order term preserving 0(3) symmetry, with 
the symmetry breaking scale v of the more fundamental 
theory. Varying the action with respect to the fields, 
we obtain following equations for scalar fields <I> m and 
neutrinos \I/ n : 

dV 
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(h Vn+ 2 f^ n <i> m = 0, (3) 



(4) 



where the 7 a -matrices satisfy the Clifford algebra in a 
locally flat inertial coordinate; 



with 



{ 7 Q ,7 b } = -2?7 ab 
if h = Diag(-1 , 1, 1, I), 



and the covariant derivative 

V a = e' i (^+r M ) 



(5) 
(6) 
(7) 



is constructed from the vierbein a and spin connection 
T M , which we give explicitly in Eqs. (25)-(31). 

Using the standard definition of the energy- momentum 
tensor 
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The energy-momentum tensor allows us to construct the 
Einstein equation, 


T l 

x ran 

Others = 


Xia 2 r . 

2r ^ 

0. 



(<JV, r + ( 5' ! lr a' ! ), 



^ + (a" 2 -!)(*„ 



(20) 



where G Ml/ is the Einstein tensor 



(10) 



(ii) 



with Ricci tensor and k — 8irG. 

When wc consider the following potential of scalar 
fields 



(12) 



with a constant 77, the global 0(3) symmetry is spon- 
taneously broken to U(l). The scalar field configuration 
describing a global monopole is known to be given by the 
hedgehog ansatz, 



$™ = F(r) — , 
r 



with a real function F(r) and r = (rEm^) 1 / 2 = (x 2 + 

y 2 + z 2 ) 1/2 , 

Assume the line element of the spherically symmetric, 
static spacetime as 



dr 2 
ot 2 (r) 



ds 2 = -S 2 (r)a 2 (r)dt 2 + + r 2 d6 2 + r 2 sin 2 0d(j) 2 



(14) 



In the Cartesian coordinates, 

x» = (t, x m ) = (t, x, y, z), for m = 1, 2, 3, (15) 
the line element can be written by 

ds 2 = g^dx^dx" ', (16) 
with the metric tensor 

= -<5 2 a 2 ^^+[^+(^-1)^^]^^", (17) 

and its inverse 

<T = -^^^+[<W+(a 2 -l)^]^ ro <T n . (18) 

From the standard definition of Christoffel symbols 

T af3 = ^"(g™, + 9v(i, a ~ ga/3, v), (19) 



In the metric given in Eq. (14), scalar field equation 
reads 



a 2 F H F ( 



(5 2 a 2 

>r i o ' r 



<5 2 



+ 2- 



cr 



4 X 2F 
+ -)- — 



w 2o„— 

aF v 

where F = 2£ F = 

We rewrite the line element as 



ds 2 = r? afc e a e fc 



with inverse of 77 , r/af, = Diag{—\, 1, 1, 1) and 

e° = a8 dt, e 1 — — , e 2 = r <f0, e 3 = r sinO dep. 
a 



(13) Putting 



& — & 11 dx . 



we find the vierbein given by 
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and their inverses 



e b = 
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The vierbein satisfy the relations 

e v e a — t) u 



(21) 



(22) 



(23) 



(24) 



(25) 



(26) 
(27) 

(28) 
(29) 



The spin connection which appeared in Eq. (7) is defined 
as 



(30) 
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with 



(31) 



We calculate the spin connections from the above Eq. 
(20) and Eq. (30), and obtain 



1 



^ 1 1 2 a cos9cos(j) 1 2 3 cos9sinq 

x 2 ' 7 r 2 7 7 rsin9 

1 o i a sincb 

+ ;t7V ' 

2 r 



(32) 



1 i o a cosOsincj) 1 7 „ cos 

r, = 2 7 7 — ; — + 2 7 7 ^r 
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2 r 

1 , ? asin9 

2 r 

which satisfy the Cartan structure equation, 

[7",IV]=7V ( 33 ) 



III. SPHERICAL SYMMETRY AND j = | 
FERMIONS 



We can show that neutrino field equations given by 
*7 a e M a(cV + r jU )* n - m„*„ + ^$ 2 *„ = (34) 
become 

n i „ i ,!,„ COS0 , 

7 + 7 *aDr*n + 7 - 09 + T^-5)*« 

ao r Ismu 

+7 3 ^^A*n " K - ^$ 2 )*„ = , (35) 
r smo v 

in the spherically symmetric and static metric in Eq. 
(14). Equivalently, we can write the last equation as 



Taking the representation of 7-matrices to be direct 
products of independnet Pauli matrices cr,-, Pj{p® o rep- 
resentation) M, 

7° = P2, 7 1 = ipi, 7 2 = -ipa<T3, 7 3 = -ip&i, (41) 
we can represent fc in 2-dimensional cr-space 



k = cr 3 — ^=d e Vsin9 + o x , - cL , 
WsinO i sinti 

and we can solve the eigenvalue equation, 

fe)r(M) = W(M), 

with 2-dimensional spinor spherical harmonics 



(42) 



(43) 



y, 



k \y>' 



(i + m)! 



1 1/2 (tanf)^/ 2 



(j - m)! 
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sin 



2JQ 



sin m 9 
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\dcos9J 2i(j-§)! 



^-(tan-) 2 Wfc , (44) 



where 



v/2|fe| 



(45) 



with k = ±(j + 5) an d j = Z ± |. In the representation 
(41) for 7-matrices, the Hamiltonian is expressed in terms 
of 2-dimensional p-matrices 



H = a5[p 3 -D r + Pl -k + p 2 (m n - ^$ 2 )] (46) 



The eigenfunctions of both H and k can be written as 
direct products of the p-space spinor 1/% and cr-space 
spinor y™, 

* n (t,r ) ^^) = e- iB " i £ tf m) «X"(M), (47) 

-j<m<j 



where the Hamiltonian is defined as 



with 



(36) 



H = c^hV-A. + ~f 1 -k + 7 °(m„ - ^$ 2 )1 , (37) 
« ir v 



IIS 2 a 2 
Dr=d r + - + -(-f + ^f). (38) 



and the neutrino equations idt^ n = H^ n then read 

r 



-^< fem ) = P3-A^ m) + [pi-+P2(m„-^$ 2 )] ^ m) - 



v 



Moreover, if we put 



(48) 
(49) 



Here we have defined the operator k as|7j H 

fc = l 7°7 1 (7 2 ^=^V^ + 7 3 ^^), (39) 

which commutes with the Hamiltonian operator H 

[k, H] = 0. (40) 



D r ^ km ^ = -(a 2 S 2 )^ d r R { n km \ 



(50) 



then we get 



aS 



rW = /03 ^a r iii &m ) + [p 1 -+p 2 (m n -^$ 2 )]^ m ). 

(51) 
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Setting 



/ (km) 
p(fem) _ I ¥>n 
rt n — (km) 

\ Xn 



(52) 



we obtain the coupled equations, 

{ ^ +iadr)ip (km) = ^_ z(mn _^|2 )]x ( fcm)(53) 

iad r ) X ^ = ^ +i (m n - 9 ^ 2 ]^, 
ad r v 

asymptotic solutions to which are given in the following 
section. 

Einstein's equation G^ u — nT^ v can be rewriten as 



(54) 



Substituting Eq. (9) into the Einstein's equation in Eq. 
(54) gives 

Rp V = K[d^ m d v $ m + g^V + l -g^ (55) 



with 



£ = e a c (*„7 c 0a* n - 0a*n7 c *»)- (56) 



Elements of Ricci tensor calculated from metric coeffi- 
cients in Eq. (14) and the ansatz for $ m in Eq. (13) are 
substituted into Eq. (55), to give the set of equations 

J_,, 2 2 ^ (S 2 a 2 ), r S 2 , r l (5 2 a 2 ), r 
2S 2{ a hrr 45 2 S 2 r S 2 

K[-V - ]p + ^(^nl dt^n ~ at* n 7°*n)],(57) 



— (& n 2 \ I ( S2a2 )>r$ 2 ,r 1 2 

2S 2[ hrr AS 2 S 2 r ,r 
n[a 2 (F, r ) 2 + V 



(58) 



1 1 m 2 t) 2 

— (1 - a z ) - -a". 

«[^ + ^ + 7 S 



2r <5 2 



2r 



(^„7 2 9 *„-a e *„ 7 2 *„)], (59) 



-(l-a 2 )--a 2 , 



1 o a 2 <5 2 , r 
r ' 2r S 2 



F 2 i 



2r sinO 



(* n7 3 9 *„-9 *„7 3 *„)]. (60) 



The first in above equations is the (t, t)-component of 
Eq. (55), the second (r, r)-component, the third (9,9)- 
component, and the last (<fr, </>)-component, respectively. 

(9, #)-component of the Einstein equation is same as 
(</>, 0)-component of Einstein equation if 

*n7 2 ^*„-a e *„ 7 2 *n = — e (^nl Z d^ n ~d^ nl H n ). 

(61) 

Using p <g> a representation for 7-matrices and the math- 
ematical identities; 



a 3 d e yr = (±<J 3 m cot9 ± vi^gW 



±a 3 e^V(j ±™ + 1)0' T m)y^ ±1 (62) 
this condition becomes 

J2 ^ (fcm 'Vi^ fcm) x 



m,m' 



(63) 



mil 



±^ ro V 3 e^ V0' ± m + l)(j T m)y f 
Tyi m ' ± V 3 e ±4 V(j ± m' + l)(j T m')K 

= e ^ (fem 'Vi^ fcm) (™+m')yr ffl -^y 

m.m' 

The condition (63) can be satisfied if 

m = to', 
to + to' = 2to = ±1, 
1 

3 - 

since 

±i±l 

Moreover, explicit calculations show us that 



(64) 



(65) 



(66) 



1 

47T' 

These make Eqs. (57)-(60) 9- and ^-independent, which 
is consistent with spherically symmetric metric in Eq. 
(14), when we consider only (j = |, m = \) or (j = 
^, m = — 2) fermions. 

For the case (j = ^, m = |) or (j = i, to = — |) with 
k = ±1, neutrino wave function is 



iE " t #' ±i) (r)}'^(^). (67) 



Let ^„(r) = tf' ^(r) or ^ fe ' *>(r) and 



(68) 
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In this case the term E in Eq. (56) reads 



with 



% 2 2 



r 2 a 2 5 2 r 2 5 r 2 aS ' 



2E n 2 ± 2i£„ 



9i = ^r 2 ad{iplp 3 d r Tp n - d r iplp 3 ip n ) 



(69) 



(70) 



(71) 



Analysing Eqs. (74)-(77) with the assumption that g,'s 
in Eqs. (70)-(73) are O(r ), we can get an asymptotic 
solution for another component of the metric in Eq. (14) 
as 



A 2 ln(— ), 
K 2M h 



(80) 



with a constant A. In this limit with the help of Eq. (68) 
the neutrino equation (53) becomes 



d r 2 ip n ~ ~q 2 {r)ip n , 



ith 



S 2 a n 4 



(81) 



(82) 



9fc = ^-raS-^-lplp^n = -^—-^(XnV + VnXn), (72) 



2tt |fc| 



2nr |fc| 



— rn — SjLn 2 



and m 

In the region 



1 i 
92 = ■^r 2 a5^] l p 2 ip n = J~(Xn<P ~ <PnXn), (73) 

Using the last equations Einstein's equations are given 
by 



— (S 2 a 2 ) ' ''" " 
2S 2[ hr 



(S 2 a 2 ), r S 2 , r | l (5 2 a 2 )„ 
4<5 2 8 2 r 



S 2 



r tt j 1/90 9l . 9fc si 

4r 2 5 2 a 2 <5 a<r 



(74) 



a 2 5 2 



r <5 2 



'"- K [a 2 (F, r ) 2 + -L(-|^ + ^)], (75) 



2r 2 M 2 a 2 



a 2 S 2 a 2 4 2F 

a 2 F +— F ( LI 1 2 ) 

2 or r r A 

dV 2g y q 2 



dF v r 2 Sa 



F = 0. 



(77) 



IV. ASYMPTOTIC SOLUTIONS FOR LARGE r 

As the global monopole solution |Q, one component 
a(r) of the metric in Eq. (14), far away from the galactic 
core, can be asymptotically taken as 



2 2 

a ~ a n = c ■ 



with a constant c when 



2M 



au 



(78) 



U(r/ 2 )~0, —(r, 2 )~0. (79) 



0<M— )< 



(83) 



~ A 2 a 2 m 2 . n 
we have an asymptotic solution to Eq. (81), 

x „ ~ ^[A n (^-a o(Z (r))e+^ (84) 

+ B n A+a q(r))e-^ r >}. 

With above asymptotic solutions to Dirac's equation for 
neutrinos, the qi in Eqs. (70)-(73) are given by 



9o 



— [B 2 n(l 



1 



-( 



^Ts 2 



■«oV(r))) 



1 /, 2 ^ 1 2 a 2 5 2 , r F 2 1 g 9l 

(76) 

and scalar field equation 



+ A n B n (l + m 2 . n )cos(2qr)), 



_ 2£, n D 2 9 2 

9l — -D n-p, 

7T 



9a- 

92 



A n B n -^-sin(2qr), 



7rr|fc| 



Trm .nC(oS 



[A n B n cos(2qr) + B n 2 ], 



with real constants A n and i?„ such that \A n \ = \B n \. 
The QiS are O(r ) or less and so are consistent with the 
assumption of Eq. (80). For r > r with In(-^fr) = 

A^ajln^ ~ > ?( r ) m Eq- (82) is pure imaginary and then 
wave functions of neutrinos become to be multiplied by 
exponentially decaying factors. Therefore we get neu- 
trino clustering with the radius r . Taking r ~ 10 23 cm 
which corresponds to the radius of our Galaxy and as- 
suming that there exists a supermassive black hole at 
the center of Galaxy with a mass Msbh — 3 x 10 6 Mq 
fll2[ , we can estimate the energy of neutrinos at E n ~ 

— • e 2 

where we have used the relations, A ^ a — — 

Mafc) and A 2 a 2 ~ 0(1). 
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V. ASYMPTOTIC SOLUTIONS FOR SMALL r 

Next let us study the small r behaviors of compo- 
nents of the metric in Eq. (14), a and 6. For small r 
near the Galactic core we adopt the Thomas-Fermi(TF) 
approximation for the energy density and the pressure 
of neutrinos, as in the case of fermion stars M ||l(j| [jTT| . 
From the local conservation law, = (^7 M \l/) ; ^ = 
^7=^^»(V — 9^1^)1 we have the conserved, total num- 
ber of neutrinos, 

N = J (Px^g^rf-® = J dPxy/^g&Ve*,). (85) 

The number density of neutrinos in the TF approxima- 
tion is given by 



d 3 q n q 



3tt 2 



(86) 



(2tt) 3 

where the Fermi distribution n q — 9(qF — q) with the 
fermi momentum qF- From Eq. (86) and the time- 
independent Dirac's equation, idt^ n — E n ^ n , we can 
put 

<^(*«7°5 t * n - 5 t W n7 °*„)) TF 

d s qn q E n {q)e\=p{r), (87) 



(2tt) 3 

where p is the energy density of neutrinos and e*o = ip- 
The stress tensor of neutrinos has the following diagonal 
elements in the vierbein basis: 



and 



1 



Pi, 
P2, 



(88) 



(^r— 9 (^j 3 d^ n ~d^ nl 3 * n )) TF =p 3 , (89) 

respectively. 

We assume, in the spherically symmetric spacetime we 
consider, that p\ — P2 — P3 = p(r). From Dirac's equa- 
tion (36) and its Hermitian conjugate, we have the rela- 
tion; 

^-(¥ n 7 ft*„-ft¥„7°*„) 
da 



(90) 



;(#n7 3 0**n - d^ n -y 3 ^ n ) = 2m eff (r) * n * r 



with m e g(r) — m n — ^-F 2 (r). Since the left hand side 
of above equation is same as i£ in Eq. (56), we have the 
following relation in the TF approximation. 



where we have defined j^-Q{r) = (^ n ^n) TF - Thus in 
the TF approximation the Einstein equations (57)- (60) 
read 



J_ r , 2 2^ (S 2 a 2 ), r S 2 , r l (S 2 a 2 ), r 
2S 2{ a hrr A5 2 S 2 r 5 2 

= K[-V + -p+-p], 



(92) 



(S 2 a 2 ), r S 2 , r 1 2 

45 2 S 2 r a " 
1 1 



= K [a 2 (F, r Y + V + -p--p} 



(93) 



: {l-a 2 ) 



1 



-a 



a 2 5 2 , r 



r" ' r 2r S 2 

, p2 „ 1 1 1 



(94) 



which are consistent with those in the case of fermion 
stars |l. The scalar field equation (21) reads 



a l F. 



. __ p ( H 

' 1 2 ^2 



a r r 



(95) 

If we assume that there exist regular solutions of the 
metric components, scalar fields, and the energy density 
and pressure of neutrinos to above equations (92)- (95) 
for small r, take series expansions as 



a 2 = a » r ™ > 52 = b n 

n>0 



n>0 



n>0 

P = 5>„r», p=J2Pnr n , (96) 

n>0 n>0 



and substitute these expansions into Eqs. (92)-(95), then 
we get the following values for the coefficients. 

«o = 1 , 01 = , a 2 = -n[^p a + ^t] 4 + ^i 2 ], 

h = , b 2 = b Q n[^p + ip + ^Vi 2 } , 

v 

Va = 0, 772 = 0, m„- — Qo = Po - 3p ■■■ 

We thus have asymptotic solutions as 
F = r;ir + C(r 3 ), 

a 2 = l- K [lp + ^+^ m 2 y + O(r 3 ), (98) 



(97) 



S 2 a 2 = {l + K [lp + l po - A„4 ]r2} + 0(r 3 )) 



where we have reparametrized the time as -y/So^ — * 
t. These solutions are consistent with the results in 



fermionic stars with a global monopole 1 1 , even if they 
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p| take F ~ for small r. The value of 771 is determined 
by the matching condition that F(r) should be continu- 
ous at r — r c , which is the upper bound for small r. This 
gives us the relation, 

Vi r c = V- (99) 

We might put r c — 2Msbh where M$bh is the mass of 
the supermassive black hole at the center of our Galaxy 

& 

VI. NEUTRINO CLUSTERING 

In Z-burst models to explain the continuation of the 
cosmic ray spectrum beyond the GZK cutoff, Weiler || 
and Fargion Q assumed that the relic neutrinos cluster 
in galaxies a few times the normal relic density, and found 
that the required flux of cosmic ray neutrinos is larger 
than previously suggested. Blanco-Pillado et al. 
suggested that the relic neutrino density in the clustering 
might be 10 12 — 10 14 times the mean relic density, tim — 
54 cm -3 . Independently in Ref. [5] a neutrino cloud was 
considered as a sphere with no diffuse boundary such that 
the relic neutrino density is given by n u (r) — nji 8(R, — 
r) + n M 0(r — R), where n R ~ 10 12 — 10 16 cm -3 and 
R ~ 10 14 — 10 20 cm. This can moderate the required 
incident flux of ultra high energy neutrinos. 

From the results in the previous sections IV and V, we 
can have a cluster of neutrinos with a diffuse boundary 
such that the density of the relic neutrinos is 

n u {r) = n c 9(r c - r) + n a {r) 6{r - r c )0{r o - r), (100) 

where n c is the constant neutrino density for small r (< 

2 

r c ~ 2Msbh) as in the section V, and n {r) = n c ^pr 
is the neutrino density for r c < r < r as in the sec- 
tion IV. Here r a ~ 10 23 cm(; the radius of our Galaxy), 
r c ~ 2Msbh — 10 12 cm(; the Schwarzschild radius of a 
supermassive black hole at the center of the Galaxy) , and 
n c >10 3 (^ff) 3 cm" 3 |§. 

VII. SUMMARY AND DISCUSSIONS 

In the most general static spacetime with spherical 
symmetry, we explicitly show that only j = 1/2 fcrmions 



can satisfy Einstein's equation, in the p ® a represen- 
tation of 7-matrices . It is also possible to show the 
fact in other representations of 7-matrices. Consider- 
ing neutrinos coupled to 0(3) scalar fields $ via <I> 2 'I"I'- 
interaction term, we had a global monopole solution jjj of 
scalar fields and asymptotic solutions of the metric com- 
ponents to Einstein's equation in Eqs. (78) and (80), for 
large r. The asymptotic solution of one metric compo- 
nent 8 2 can be given in series of more general functions as 
$ 2 {r) ^ J2n fn 2 (r) lnn (jM) witn fn(r -> 00) = constant, 
which do not change the asymptotic behaviors of other 
solutions in this limit. 

We adopted the Thomas-Fermi (TF) approximation for 
small r expansions of the energy density, pressure of neu- 
trinos and so on. Since only s-waves are considered in 
the TF approximation, small r expansions in the sec- 
tion V are consistently related to asymptotic solutions 
obtained for large r in th esection IV. More rigorous con- 
nection between two asymptotic solutions can be made 
by numerical methods. With the simple assumption, Eq. 
(100), about the radial dependence of the relic density 
of neutrinos clustered in the galactic halo, further stud- 
ies of Z-burst models shall be possible. For small r so- 
lutions in the section V to have astronomical meaning, 
the global 0(3) symmetry breaking scale n and the scalar 
self-coupling constant A in Eq. (12) should be very small, 
which might be realized in some Majoron models ]T^ |. In 
such a case the deficit angle 1 — c = nn 2 in Eq. (78) is 
negligible M. 
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